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We examine the Friedel sum rule which states that the "excess charge" due to a single 
impurity potential in a metal is equal to a sum of phase shifts for scatterings of electrons 
by the impurity. For finite volume, the "excess charge" is given by the difference between 
total numbers of levels in the Fermi sea with and without the impurity potential. However, 
a sequence of the "excess charge" for finite volume is not necessarily bounded in the infinite 
volume limit, as was pointed out by Kirsch. In order to circumvent this difficulty, we define 
"excess charge" directly for the infinite volume. The Friedel sum rule is proven to hold 
for the "excess charge" thus defined. 
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1 Introduction 

In this paper, we examine the Friedel sum rule [3, 4, 7]. Consider a metal with a single 
impurity at zero temperature. Electrons are scattered by the impurity potential, and the 
charge distribution of electrons below the Fermi level changes. "Excess charge" is defined 
to be the difference between the total numbers of levels in the Fermi sea with and without 
the impurity potential for a fixed Fermi energy Ep. Then the Friedel sum rule states that 
the "excess charge" can be expressed in terms of a sum of phase shifts due to the impurity 
potential. However, a sequence of the "excess charge" for finite volume is not bounded in 
the infinite volume limit, as was pointed out by Kirsch [6]. Although this result may not 
necessarily imply that there exists no sequence of finite boxes such that the sequence of 
the corresponding "excess charges" converges to a finite value in the infinite volume limit, 
it clearly implies that the "excess charge" cannot be well-defined in the infinite-volume 
limit from a sequence of finite boxes. In order to circumvent this difficulty, we define 
"excess charge" directly in the infinite volume. We prove that the "excess charge" thus 
defined can be expressed in terms of a sum of phase shifts of scatterings of electrons by 
the impurity potential. Thus we justify the Friedel sum rule in the infinite volume. 

In the next section, we describe the model and our results. The proofs of the main 
theorems are given in Sections 3 and 4. 



2 Model and the results 

Consider a single electron in an impurity potential V on R 3 . The Hamiltonian H is 

H = -A + V. (2.1) 

We assume that V is bounded and of compact support. Let P be the projection on energies 
smaller than the Fermi energy Ep. We denote by Hq = —A the Hamiltonian without the 
impurity potential, and denote by Q the corresponding projection on energies smaller than 
the same Fermi energy Ep. We also denote by x = £2, £3) G R 3 the coordinate of the 
electron. Let xr De the characteristic function of the ball defined by B R = {x e R 3 | \x\ < 

R} centered at x = with the radius R > 0, where |x| := \Jx\ + x\ + x\. We define the 
excess charge by 

Z := HmTr Xr{P ~ Q)XR- (2.2) 

KjOO 

Theorem 2.1 The limit, lim^^Tr xr(P — Q)xr> exists for all the Fermi energy Ep > 0. 
The proof will be given in Section 3. 

The Friedel sum rule [4, 7] gives a relation between excess charge Z and the number 
of bound states and phase shifts for a spherical potential V. The precise statement is as 
follows: 

Theorem 2.2 Assume that the impurity potential V is spherical, i.e., V is a function of 
the single variable \x\. Let Nb(Ep) be the number of the bound states of the Hamiltonian 
H with the energy smaller than the Fermi energy Ep > 0, and let rjz{k) be the phase shifts 
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of the scattering states with angular momentum £ and with wavenumber k. We choose the 
phase shift r)g(k) to be continuous in all k with lim^oo f]i{k) = for all £. Then 

1 oo 

Urn Tr X r(P - Q)xr = N b (E F ) + - £(2£ + 1) [ m (k F ) - m (0)\ , (2.3) 

where kF is the Fermi wavenumber, i.e., k? = y/Ep. 

The proof will be given in Section 4. 

Remark: I. Since the Friedel formula (2.3) states that the difference between the two 
bulk quantities is related to the boundary quantities, the formula can be geometrically 
interpreted as an analogue of the Gauss or Stokes theorem . 

2. In one and two dimensions, one can easily obtain a similar formula to the Friedel sum 
rule (2.3) for the phase shifts. However, our method cannot be applied to the systems in 
higher dimensions d > 3 because of the singularity of the Green function at short distance. 
See eq. (3.8) and the corresponding argument in the next Section 3 for details. 

3. For the present spherical potential V, one can prove the absence of bound states with 
the binding energy E = and with the angular momentum £ = 0. Further, the Levinson 
theorem [8] implies that 1 

r?0 (0)=0 for £ = 0, and (2£ + 1)^(0) = nN e for i > 1, (2.4) 

where N# is the number of bound states with angular momentum £ > 1, including those 
with zero binding energy. Substituting these into the above result (2.3), we have the 
well-known form of the Friedel sum rule [4, 7] as 

hm Tr X r(P - Q)Xr = ~ E( 2 ^ + X MM (2.5) 

i?too TT e=Q 

for the Fermi energy E-p > 0. 

4. We should remark the relation between the excess charge Z and the "Krein spectral 
shift" . Actually the excess charge Z is equal to the spectral shift in a formal calculation. 
For the spectral shift, see recent papers of Combes-Hislop-Nakamura [1] and Hundertmark- 
Simon [5] and references therein. 

3 Proof of Theorem 2.1 

We begin with recalling the well known fact that the spectrum a(H) of the present Hamil- 
tonian H of (2.1) is absolutely continuous in the region of the positive energy. See, for 
example, the book [10]. The projections below the Fermi level E F > are, respectively, 
written 

„ 1 f dz , ^ 1 r dz , 
p = — and Q = — / — 3.1 

2mJrz-H ^ 2m Jr z - H y J 

1 See, for example, Chap. 12 of the book [9] of Newton. 
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by using the contour integral and the resolvent. From this expression, one has 



^ 2m Jr 



dzl- 1 ^-- 1 ^)- (3-2) 



,z-H z-H Q 

The difference of the two resolvents in the integrand is computed as 

1 1 1 -v- 1 



z — H z — H z — H Q z — H 

= — - — V — - — + — - — V—^—V — - — . (3.3) 

Z — Hq z — Hq z — Hq z — H z — Hq 

These observations yield 

If 11 

Tr Xr(P ~ Q)Xr = 7T- dzTT XR t-V — Xr 

2m Jr z — H z — H 

+ -L/ dz Tr Xr ^—V-^—V—^—xr- (3.4) 
2m Jr z — H z — Hz — H 

It is not difficult to evaluate this right-hand side [11] except for the contributions of the 
contour integral near the Fermi level. In the following, we treat only the contributions near 
the Fermi level. Therefore we write z = (k-p + ifi) 2 near the Fermi level Ep by introducing 
the real variable fx. Then the integral kernel of the resolvent (z — Ho)^ 1 with z = (kp + ifi) 2 
is written 

G (x; (k F + ifx) 2 ) = -— — -exp[(ik F - n)\x\], (3.5) 

4:7l\X\ 

and the contribution of the contour integral near the Fermi level is written as 

If 1 r-e r+a 

/ dz- ■ ■ = lim- / +/ dux (k F + ia)--- (3.6) 

2m J elO 7T J - a J+e 

with a positive cutoff a. 

Using the representation (3.5) of the Green function, we have 

Tr xr — — l ~ ir XR 
z- Ho z- Ho 

= Y^J d3x l d " x ' Xr(x) j-p exp [2(ifc F -n)\x- x'\] V{x'). (3.7) 

In passing, we remark the following: In a higher dimension d > 3, the Green function 
G (x; z) behaves as 

„ , , Const. 

G (x;z)~--^ (3.8) 

for the short distance \x\ ~ 0. In that case, the right-hand side of (3.7) becomes a divergent 
integral. Therefore our method cannot be applied to the higher dimensional systems. 

We return to the present three dimensional system. In this case, the contribution for x 
in finite range in the integral of the right-hand side of (3.7) gives a finite value. Therefore 
it is sufficient to evaluate the following integral: 

h — I {k F + ijj)djj, [ d 3 x- -y- exp [2(ik F — jj)\x — x'\] , (3.9) 

Je JB R \B R , \x — x'\ l 
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where we have chosen the radius Rl of the ball B R > to satisfy B R i C B R and B R/ / 2 D suppV, 
and we have taken account of the contour integral near the Fermi level with the cutoff 
constants a, e > as mentioned above. We will take the limit e [ 0. Note that 

pa { 

I (k F + ifj)dfiexp[— 2/jl\x — x'\] = 



2\x - x 



— exp[— 2fi\x — x'\](kF + ifJ>) 



+ 



2\x — x 



f a 

— / d/xexp[— 2/jl\x — x'\ 

Je 



h + ^exp[-2 £ |:r - x'\] + 0(\x - x'\~ 2 ) 



\x — X 



(3.10) 



for a large \x — x'\. Substituting this into (3.9), we have 

lim/i = lim / d 3 x kF + Z f q exp [2(ik F - e)\x - x\] + 0(1) 

eiO s10Jb r \B r , \X-X'\ 3 PL V ^ U K ' 

k F 



= I d 3 x-^—^exp[2ik F \x-x'\] + 0(l). 



(3.11) 



By taking the radial variable r = \x — x'\ centered at x', the integral about the radius r 
in the first term can be written 

roo 1 

/ dr-e 2lkpr (3.12) 
Jro r 

in the limit R j oo, where i? is a cutoff constant given by i? — ^o 71 "^ 1 with a positive 
integer n . Since both of the real and imaginary parts of the integral can be treated in the 
same way, we consider only the imaginary part. This integral converges to a finite value 

as 



f°° , sin2£; F r 
/ dr 

JRo 



= E 

oo 

= E 

m=no 



r «+i sin2A; F r 
dr 



''2 



r 2 m+i sin2/c F r / >r 2m+2 sin2/c F r 

r 

1 



h / dr- 



~ pm+1 / 1 

= E / ^ r srn 2A; F r - 

m=n "' r 2'» \ r 
/-r-2 m +i 

> / sin2/c F r— — t— — 

±^,Jro m r(r + 7rkp l /2) 



r + nk^/2 

TT/tp 1 / 2 



< oo, 



m=no 



(3.13) 



where r n = nnk F 1 /2 with the integer n > n . 

The integrand of the second term in the right-hand side of (3.4) is written 



Tr XR —!—V—^-V—^— X r 



' z — Hq z — H z — if o 
= / d * x J d * x ' \ d 3 x"x R (x) _ exp[(ife F - fi)\x - x'\] 

x V(x')G(x',x"] (k F + ifi) 2 )V(x")—^ — r exp[(ik F - n)\x" - x\], (3.14) 

Ob Ob 
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where G(x,x'; (k F + i/i) 2 ) is the integral kernel of the resolvent (z — if) -1 . Since the 
contribution of the integral for x in finite range is finite, we consider only the integral for 
x in the region Br\Bri. We choose R' large enough so that the ball B R i/ 2 includes the 
support of the potential V . Then we have 

/ d 3 x I d 3 x' I d 3 x"- — — - exp[(ik F - /i)\x - x'\] 

x V(x')G(x', x"; (hp + ifi f)V{x")—^ - exp[(iA; F - fi)\x" - x\] 



( d z x^-,e 2 ^-^ x \ (u, Viz - if)" W 
Jb r \b r/ \x\ 2 \ / 

/ d 3 x-^e 2{lkp -^ xl {u(x,-),V(z-H)- 1 Vd(x,-)) (3.15) 
JBr\b b i \x\ 



by using 



and 



\x — x'\ \x\ 



J +0(\x\- 2 ) (3.16) 



\x-x'\ = \x\ -^- + 0{\x\~ l ) (3.17) 



\x\ 

\2 



for a large \x\. Here z = (k F + i/i) , and the functions u and v are given by 

u(x') = exp[(ik F + n)x ■ x' /\x\\ (3.18) 

and 

v{x") = exp[(-ik F + n)x ■ x" /\x\], (3.19) 

and u(x, x') and v(x,x') are bounded functions. Note that ||xa(<2 — -£0 _1 Xa|| with z = 
(k F +ifi) 2 is bounded uniformly [2] in fi > for the Fermi energy E F > and for a compact 
subset A of R 3 , where xa is the characteristic function of A. From this and the method 
in (3.10), the contribution form the second integral in the right-hand side of (3.15) can be 
handled in the same way. Thus it is sufficient to evaluate the integral, 

I 2 = f a (k F + ifi)dfi f d 3 x-^ exp[2( ik F - v)\x\](u,V(z- H^Vv) , (3.20) 

Je JB R \B Rl \X\ 2 X ' 

with z = (k F + in) 2 . Integrating by parts, we have 

f 1 (i 

h = ~{k F + ifi) d 3 x—-^exp[2( y ik F -fi)\x\](u,V( y z-Hy 1 Vv) a 
Jb r \b r , 2\x\ 6 \ I £ 

+ d/2 d 3 x 7 ^exp[2(ik F -i2)\x\} — (k F + i/2)(u,V(z-H)- 1 Vv). 

(3.21) 

The first term is written 

-(Jfep + ie) [ d 3 x—^exp[2(ik F -e)\x\}(u,V(z- H)- l Vv) + 0(1) (3.22) 
Jb r \b r i \x\ 
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with z = (k F + is) 2 . Note that (u, V(z — H)~ l Vv) does not depend on the radius \x\ and 
is bounded uniformly in e > for the Fermi energy E F > 0. From these observations and 
the same argument as that for the integral of (3.11), we have that the first term in (3.22) 
converges to a finite constant as R j oo and e j 0. 

Similarly, the second term in the right-hand side of (3.21) is written 

- f 1 dp(k F + ipf [ d 3 x-^—exp[2( y ik F -p)\x\}(u,V( y z-H)' 2 Vv) + 0( y l). (3.23) 

Je JB R \B R , 2\x[ i X 1 

Therefore it is sufficient to show that | (u, V(z — H)~ 2 Vv) \ is bounded uniformly in /j, > 0. 
Using the resolvent identity, one has 

V— — V = V — —V + V^—V — — V. (3.24) 

z-Hz-H z-H Q z-H z-H z-Hqz-H 1 ; 

Moreover, one obtains 

V — — V = V — — V + V — — V — - — V. (3.25) 

z-H z-H z-H z-H z-Hqz-Hq z-H y ' 

From these two identities, it is sufficient to show boundedness of the operator, 

V 7 l H 7 l H V- (3-26) 
z — n z — n 

Using the spectral decomposition, one has 

= (2^ £_r dp l du v ~ (hi + i^ f ^ u)p2 + 

(3.27) 

for two functions, ip, ip, of compact support, where we have introduced the polar coordinate 
k = puj with p = \k\, 5 is a positive cutoff, and the function / is given by the product of 
the Fourier transforms of ip,ip. Note that 

1 19 1 , 

(3.28) 



[p 2 - (k F + ip) 2 } 2 2pdpp 2 -(k F + ip) 2 
Using this identity and integration by parts, one has 

rk F +8 I I rk F +6 I Q 

Ls dp [p 2 -(k F + i p) 2 ] 2fMp =2jU dp 7^ihTWdp- fMp + 0{1) - 

(3.29) 

The integral of the first term in the right-hand side is written 

rk-F+5 \ _ rk F +5 \ . _ _ -, 

L-s dp p 2 -(k F + l p) 2 f{kF > u) + Ls dp p 2 -(k F + l p) 2 [fM-f^l (3.30) 
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where we have written 



d 

f(p,u) = —f(p,u)p. 



Note that 



and 



1 



p 2 — (kp + ip) 2 2(kp + ip) 



P 



p — kp — ip p + kp + ip 



ip 



(3.31) 



(3.32) 



(3.33) 



p — k F — ip (p — k F ) 2 + p 2 (p — kp) 2 + p 2 

Combining this with the differentiability of the function /, we can prove that the right- 
hand side of (3.29) is bounded uniformly in p > 0. Thus the double limit R | oo and e I 
for the integral (3.20) exists. Consequently the limit, lim^oo Tr xr(P ~ Q)xr-> exists for 
the Fermi energy E F > 0. 



4 Proof of Theorem 2.2 

We give the proof, following the idea of Friedel [4]. We denote by J := —ix x V the 
angular momentum, and define a regularized excess charge as 

Z PtR :=Tie-^ X R{P-Q)XR ^r [3 > 0. (4.1) 

Since the integral kernel of e _/3j2 which is the heat kernel on the two-sphere S 2 depends 
only on the angular variables except for the identity about the radial part, the existence of 
liniR^oo Zp^R is justified in the same way as in Section 3, and we recover the excess charge 
Z of (2.2) in the limit (5 { as Z = lim^ lim^oo Zp R . 

The integral kernel of the projection P for the Hamiltonian H with the impurity 
potential V is written 



P(x, x') = Yl v n{x)v*{x') + - f kF dk ]T u kAm {x)u k/ * m {x') (4.2) 

n:E n <E F 71 J ° i,m 



in terms of the bound state v n and the scattering state Uk/, m with the angular momentum 
£ and the magnetic quantum number m. Similarly, 

Q(x,x') = - [ F d*£« M ,£WM,m V) (4-3) 

in terms of the scattering state u k J® of the Hamiltonian H . The scattering state 
is written 

^ m (x) = ^ ra (^) (4.4) 

in terms of the angular part Y™ and the radial part f k j. The radial part follows from the 
equation, 



d f 



k 2 - V(r) - 



+ 1 



/M = °- ( 4 - 5 ) 
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Multiplying this by another solution fy ^ with a different wavenumber k', and then inter- 
changing the variables k and k' in the resulting equation, and finally taking the difference 
between the two equations, one has 



d 
dr 



f , —f -f —f, 



+ (k 2 - k ,2 )f k , : ef k:e = 0. 



(4.6) 



Integrating over r from to R, one obtains 

d „ „ „ s „ , „, d 



fk>,e(R)-^fkA R ) - fkAR)^fk'AR) = (k' 2 - k 2 ) J q drf k , t e(r)f k Ar), (4.7) 

where we have used the fact /^(O) = 0- Further, by dividing by k' — k and taking the 
limit k' — > k in the both sides, one has 



I R -1 



R 



dr \fkA r ) 



(4.f 



We will justify differentiability of the function with respect to k later. The function 
fkA r ) f° r a large r outside the support of V(r) = is written as 



fkA r ) = ^2 [krjf(kr) cosr]g(k) — krng(kr) smrjg(k)] 



(4.9) 



in terms of the Bessel functions jg and ng. The asymptotic forms of the Bessel functions 
for £ > 1 are given by 



m = - p 



sin(p - n£/2) + cos(p - n£/2) 



and 



(£ + IV 1 

cos(p - ir£/2) - )„ ,'. — sin(p - ir£/2) 



(e-i)\2 P 

for a large p. This yields the suitable asymptotic form, 



(4.10) 



(4.11) 



hAr) ~ sin(A;r + rp(/c) - vr£/2) + |±1|-^- cos(A;r + Tfc(fc) - n£/2), (4.12) 

for £ > 1. The case with £ = is much simpler. Substituting the expression (4.9) of f k ,g 
into the left-hand side of (4.8), one has 



jR 

1 k,e 



dr)g(k) 



+ Rcos 2 7] e (k) 



(j'g(kR)) 2 -jg(kR)~j';(kR) 



dk 

+ Rsin 2 r]g(k) [(n' e (kR) f - n e (kr)n' e \kR) 

- Rsinr]g(k) cosrjg(k) [2j' t (kR)n' t (kR) -jg{kR)n"{kR) - jg (kr)hg(kr) 

- i cos 2 T]g(k)jg(kR)j'g(kR) - ^ sin 2 r] t (k)hg{kR)h'g{kR) 
+ — svn.r]g{k) cosr)g(k) je(kR)n' e (kR) + n e (kR)j' £ (kR) 



(4.13) 
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where we have written ji(p) = pje{p) and n^{p) = pni(p), and used the Wronskian, 

je(p)n' e (p)-n e (p)~fe(p) = l- ( 4 - 14 ) 

Here /' is the derivative of the function /. Taking the difference between the two quantities 
lj? e with and without the impurity potential V(r), one obtains 



jf<fr{|/ M (r) 



lS(r) 



driejk) 
dk 



+ g e (k,R), 



(4.15) 



where 

9e(k,R) 



-Rsin 2 r]i(k) 



(j' e (kR)) 2 -j e (kR)%(kR) 



+ Rsm 2 r} e (k) [{h' t (kR)f - n t (kr)n'l(kR)] 

- R sin r] £ {k) cos r] e {k) \2j' e (kR)n' e (kR) - j e (kR)n"(kR) - ~j" \kr)n E {kr) 



+ 



sin 2 ri e (k) \j e (kR)j' e (kR) - n e (kR)n' e (kR) 



+ - sin r] e (k) cos rje(k) je(kR)n' e (kR) + ne(kR)j' e (kR) 

rC 



(4.16) 



As is well known, the function fk,e(r) is continuous 2 with respect to r > and to k > 0, 
and the function ge(k,R) is also continuous with respect to k > 0. Thus the existence of 
the derivative of rje(k) in the right-hand side of (4.15) is justified for k > 0. Clearly this 
implies that the function is also differentiable with respect to k > for a large r from 
the expression (4.9). We decompose the function g e (k,R) into two parts as 



ge(k,R) = ge(k,R) - ^smr) e (k) cos(2kR + r) e (k) - ir£). 



Note that the phase shift rje(k) satisfies 3 

| sin^(/c)| < Const, 
and that the Bessel functions behave as 



for £ = 0; 
i*- 1 for i > 1, 



j e (p) ~ ap e+1 + bp 



J+3 



and 



hg(p) ~ cp 1 + dp 



-1+2 



(4.17) 



(4.18) 



(4.19) 



(4.20) 



for a small p, where a, b, c, d are the real constants. Using these properties, one can prove 
that the function ge(k, R) is bounded for any k G [0, fcp]. Further, by using the asymptotics 
(4.10) and (4.11) for a large p, one has 



g e (k, R) — > as i? — > oo for a fixed 7^ 0. 



(4.21) 



2 The continuity of the function /^^ can be proved by using the method of integral equation. See, for 
example, Chap. 12 of the book [9] of Newton. 
3 See, for example, Chap. 12 of the book [9]. 
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From these observations, the dominated convergence theorem and (4.15), we obtain 

i&r dr r <ft { |/w(r)|a -i / swr} 

= Ve(k F ) ~ %(0) - lim dk ^J^ cos(2kR + rjAk) - ir£), (4.22) 

R^co Jo k 

where the phase shift with k = is defined as rje(0) := lim fc | rj e (k). The oscillating integral 
in the right-hand side vanishes in the limit R f oo by the Riemann-Lebesgue theorem and 
the boundedness of sin rje(k)/k which holds because of the bound (4.18). Combining this 
result with (4.1)-(4.4), one has 

Jim ZpR = £ e-MW + ± £(2£ + IK [ V e(k F ) - Ve (0)\ . (4.23) 

bound states £=0 

We define the sum about £ in the right-hand side of (2.3) by the limit j3 { of the second 
sum in the right-hand side of (4.23). Then the desired result (2.3) is obtained in the limit 

/no. 
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